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Abstract 



We study nuclear effects in the structure function F3 which describes the 
parity violating part of the charged-current neuitrino nucleon deep inelastic 
scattering. Starting from a covariant approach we derive a factorized expres- 
sion for the nuclear structure function in terms of nuclear spectral function 
and off-shell nucleon structure functions valid for arbitrary momentum trans- 
fer Q and in the limit of weak nuclear binding, i.e. when a nucleus can be 
treated as a non-relativistic system. We develop a systematic expansion of 
nuclear structure functions in terms of Q~ 2 series caused by nuclear effects 
("nuclear twist" series). Basing on this expansion we calculate nuclear cor- 
rections to the Gross-Llewellyn-Smith sum rule as well as to higher moments 
of F3. We show that corrections to the GLS sum rule due to nuclear effects 
cancel out in the Bjorken limit and calculate the corresponding Q~ 2 correc- 
tion. Special attention is payed to the discussion of the off-shell effects in 
the structure functions. A sizable impact of these effects both on Q 2 - and 
x-dependence of nuclear structure functions is found. 
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I. INTRODUCTION 



Experiments on deep-inelastic scattering (DIS) of charged leptons and neutrino remains 
an important source of information about nucleon and nuclear structure. Let us note in this 
respect that, though the primary aim of these experiments is to explore the structure of 
hadrons, DIS data is collected usually for nuclear targets. Since the discovery of nontrivial 
nuclear effects in DIS with muon and electron beams (the EMC effect), the studies of nuclear 
structure functions is an important topic of itself (for a recent reviews of experimental 
situation and theoretical approaches to nuclear effects in DIS see 

In the literature nuclear effects have been discussed in detail for the spin independent 
structure functions as well as for spin structure functions gi$ which are measured in the 
charged lepton DIS. However up to now only a little attention was payed to nuclear effects 
in neutrino DIS [0-0] an d especially to the structure function F3 which does not have an 
analog in the DIS with the charged leptons. Perhaps such a lack of interest can be partially 
explained by the lack of precise neutrino data for different targets, that would allow to 
separate nuclear effects. We recall also that in the parton model F 3 is related to the valence 
quark distribution and at large x > 0.3, where the sea part is small, xF 3 is practically equal 
to F 2 . It is expected therefore that nuclear modifications of xF 3 is similar in this region to 
those of F 2 , which has been studied in more details. It should be stressed however that this 
similarity between F 2 and xF 3 is only approximately valid in a certain kinematical region 
and fails at small x, where the sea part of the parton distributions is important, as well as 
when one studies effects due to finite Q 2 . We note also that a precise neutrino data is now 
available for the iron target || (for a review on neutrino data see Ref. [|J), and in order to 
perform a QCD analysis || of the data one should have a reliable method to handle nuclear 
effects in the neutrino DIS at finite Q 2 . 

These arguments motivate us to study the neutrino and anti-neutrino DIS from nuclear 
targets, that we do in the present paper. Our primary aim here is to study nuclear effects 
in the structure function F% within approach applicable at finite Q 2 and develop a reliable 
method to extract F3 of the isolated proton and neutron from nuclear data. For this purpose 
we re-examine and extend to finite Q 2 the impulse approximation formalism for nuclear 
structure functions, developed previously for large Q 2 for the structure functions F lj2 jTT| ) [12 
as well as for the spin structure functions gi% 



The paper is organized as follows. In Secs.p|,p| we present a general covariant framework 
in which nuclear structure functions are expressed in terms of nucleon propagator in the 
nucleus and the off-shell nucleon hadronic tensor. Sec. [lll B| describes how the relativistic 
expressions for nuclear structure functions can be reduced in the limit of weak nuclear 
binding by making a non-relativistic expansion of nuclear matrix elements. In Sec. [lll O we 
analyze the resulting generalized convolution equation for the nuclear structure function 
and develop its systematic expansion in nuclear Q~ 2 series. In SeafTV] we discuss off-shell 
effects in the nucleon structure function. Nuclear corrections to the Gross-Llewellyn-Smith 
sum rule as well as to higher moments of F 3 are calculated in Sec.|V[ Sec.[VT] deals with an 
approximate expression for the nuclear structure function which arises if one treats the effects 
of Fermi-motion and nuclear binding perturbatively. We summarize and make concluding 
remarks in SeclVll. 
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II. HADRONIC TENSOR AND STRUCTURE FUNCTIONS 



We consider the neutrino (anti-neutrino) charged-current interaction. Inclusive inelastic 
scattering from hadrons is described by the hadronic tensor 

W%{P, q,S) = i- J d'z e *" (P, S | [j£(z) t J?(0)] | P, S) , (1) 

where is the charged-current corresponding to the exchange of W /± -boson, P and q 
are the four-momenta of the target and virtual IV-boson, respectively, and the vector S is 
the target polarization (spin) vector. In what follows we discuss scattering of unpolarized 
neutrino beams from unpolarized targets and will drop an explicit polarization dependence 
of the hadronic tensor. 

It is well known (see, e.g., [HJ) that unpolarized neutrino nucleon deep- inelastic scat- 
tering is described in general by five structure functions. Two of those are similar to the 
structure functions F\ and F 2 in charged lepton DIS. The structure function F 3 arises due to 
vector-axial current interference in Eq.(^) and determines the v-v cross-section asymmetry. 
The left two structure functions, F4 and F$, are due to nonconservation of the axial current. 
These structure functions are usually not discussed in the DIS regime, because their contri- 
bution to the cross-section is suppressed by a tiny ratio of the lepton mass to the neutrino 
beam energy. Apart of this argument, F 4 and F$ have some "internal smallness" due to 



PCAC 0. 



In the present paper we focus on the discussion of the vector-axial current interference 
term which contribution to the hadronic tensor ([!]) reads: [] 

^(P,g)=u^PV F f p ^ , (2) 

where the four-momentum transfer squared Q 2 = —q 2 and the Bjorken variable x = Q 2 jlP-q 
are the standard DIS variables. 

In our covariant analysis of neutrino-nucleus scattering it will be useful to work with the 
off-shell nucleon tensor W^, which is defined through the imaginary part of the forward 
VF-boson scattering amplitude from an off-shell nucleon. In terms of W^, the hadronic 
tensor of a "physical" on-mass-shell nucleon averaged over nucleon polarizations is: 

W£(p, q) = ^Tr [ty + M)W, u {p, q)] , (3) 

were p is the nucleon momentum, M is the nucleon mass and p 2 = M 2 . 

In off-shell region the Lorentz tensor structure of is more complicated than the 
corresponding structure for the on-shell nucleon. In order to establish the tensor structure of 
W^j, we expand the latter in terms of a complete set of Dirac matrices, {/, 7", a a/3 , 7^75, 75}. 
The various coefficients in this expansion must be constructed from the vectors p and q, and 



1 We use the following normalization of the states throughout of the paper (p\p') = (27r) 3 2po<5(p — 
p'). In this normalization the hadronic tensor ([!]) is dimensionless. 
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from the tensors g a p and e^vap- We consider here the anti-symmetric part of the tensor 
Wfj, u and separate only those terms which are even under time-reversal and odd under 
parity transformations, since only such terms can contribute to F 3 . If we keep only current 
conserving terms then the most general form of the antisymmetric tensor W M „(p, q) which 
obeys to these requirements can be written as follows,^ 

WV„(p, q) = i E„ ra/ i 9° S"(p, <z)/(4p • q), (4) 

where the coefficient functions p l > are scalar, dimensionless and real functions of q 2 , p ■ q 
and p 2 . 

We find therefore that an off-shell nucleon is characterized by four independent structure 
functions which correspond to T-even and P-odd terms in the hadronic tensor. On the mass- 
shell p 2 = M 2 , because of the Dirac equation, only one out of the four structure functions 
in Eq.(|5|) is independent, and Eq.(|4]) reduces to Eq.(H). Substituting Eqs. (|4],|5|) to Eq.(|) we 
find: 

F g = / (0) +/ ( P ) +/ (9) +/ (7) j 

with the functions /W evaluated here at p 2 = M 2 . 

We conclude this section by noting that the nucleon structure function F3 splits into four 
structure functions /W in off-shell region. We discuss the concequences of this observation 
for nuclear scattering in the next section. 



III. NUCLEAR STRUCTURE FUNCTIONS 

Discussion of nuclear effects in DIS is usually done separately for small and large Bjorken 
x. The corresponding scale for this separation comes from the comparision of a characteristic 
time for DIS, 1/Mx (see, e.g., Ref. [0), and an average distance between bound nucleons 



which is about 1 Fm for heavy nuclei. At large x > 0.2 the characteristic DIS time is smaller 
than average internucleon distance. This observation motivates the use of the impulse ap- 
proximation in this region, i.e. the nuclear Compton amplitude is approximated by the 
uncoherent sum of the scattering amplitudes from bound nucleons, Fig.|I|. Major nuclear 
effects found in this region are due to Fermi- motion fl9" Hl3|l nuclear binding |ToHl3| , and off- 
shell modification of nucleon structure functions [Fj-inj] (for more discussion and references 
see review paper 0). 

At small x DIS proceeds mainly via diffractive mechanism, i.e. virtual photon/boson 
first converts into a qq fluctuation which subsequently interacts with the target. Since an 
average life time of the qq fluctuation at small x is large as compared with an average 



2 The symmetric part of the tensor W^, which is relevant for F% and F2 structure functions, have 
been discussed in detail in Ref. [14]. The T-even and P-even terms in the antisymmetric part of 
Wui/, which determine the spin structure functions g\ and 52 > have been constructed in Ref. [13]. 
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distance between bound nucleons, coherent multiple scattering effects are important. It is 
known that for the structure function F 2 this mechanism leads to nuclear shadowing effect 
(for a review of models for nuclear shadowing see, e.g., However a little is known about 
shadowing effect in the structure function F 3 ||. 

Other possible nuclear effects which go beyond the impulse approximation are corrections 
due to meson exchange currents (MEC) and final state interactions (FSI). 

While MEC give a sizable correction to the structure function F 2 at small x < 0.2 (see, 



e.g., Refs. [0,|TT|] and a review 0), one may argue that for F 3 the corresponding correction is 
less important. Indeed the structure function F 3 in neutrino scattering in the quark-parton 
model reads F£(x) = d(x) — u(x) + s(x) — c(x). For light mesons, such as ir- and p-mesons, 
the contribution from u- and ci-quarks cancels out because of the isospin symmetry. F3 
therefore vanishes for light mesons if one neglects a contribution due to s — c (similarly 
F 3 vanishes for anti-neutrino scattering up to s — c contribution). By this reason one can 
neglect MEC correction to nuclear F 3 . [] 

It is commonly believed that the FSI in DIS is a higher twist effect, i.e. appears at least 
in Q~ 2 order. It has been indeed found within framework of the collinear expansion fll? 



that the Feynman diagrams describing the FSI of the struck quark with the rest of the target 
are suppressed at high Q 2 . Similar conclusion is drawn within Bethe-Solpeter formalism in 
Ref. ||18|| . Cancellation of the FSI in the Bjorken limit is also found in the studies of inelastic 
scattering in non-relativistic [|T5| , p~9|] as well as in relativistic |2(J two-body models even with a 
"confining interaction" such as oscillator potential. These observations might be considered 
as an argument in favor of the "common wisdom" that cancellation of FSI effects in total 
cross-sections is due to unitarity and closure of the final states. However the situation is less 
clear in the case of DIS from nuclei, where the struck quark may scatter softly from bound 
nucleons while before hadronization. This effect may lead to a finite FSI correction to the 
nuclear structure functions even in the Bjorken limit Phenomenological arguments in 
favor of the suppression of the FSI in nuclear DIS are recently given in Ref. [ 22| . Finally we 
comment that the systematic treatment of the FSI in DIS from nuclei is still lacking and 
remains a challenging problem, which is however out of the scope of the present paper. 



A. Impulse Approximation 

Basing on the discussion given above we consider here the nuclear Compton amplitude 
in the impulse approximation (IA), Fig.[IJ. Our aim here is to discuss the IA in terms of a 
covariant framework valid at any momentum transfer and eventually develop a systematic 
expansion scheme in terms of Q~ 2 for the corrections due to nuclear effects. 

We take the imaginary part of the nuclear Compton amplitude in IA and find that the 
hadronic tensor of the nucleus with momentum Pa can be written as follows: 

W*{P A , q) = £ / \M Tr [^(P; Pa) VV^Cp, q)} , (7) 

r=p,n 



3 We comment however that MEC may cause a finite contribution to nuclear F3 due to s-quark 
content of the mesons. Special care has to be taken of the strange mesons contribution. 
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where [dp] = d 4 p/ (27r) 4 and W^(p, q) is the hadronic tensor of the off-mass-shell proton or 
neutron, given by Eqs.(f|) and (Q). 

The function A T (p; Pa) is the proton (r = p) or the neutron (r = n) propagator in the 
nucleus which is defined as follows :Q 

K P {p\Pa) = J dhdte 1 ^-*-* (P A \N;(0)Nl(r,t)\P A ) . (8) 

Here we describe nucleons by the field operator N^(r,t), the averaging is taken over the 
nuclear state with momentum P A . 

The expression for the nuclear tensor in Eq.(|7]) is covariant and can be evaluated in any 
reference frame. It will be convenient, however, to work in the target rest frame, in which 
the target four-momentum is Pa = (Ma',0), and the momentum transfer to the nucleus, 
q = (go; 0j_, — |q|), defines the z-axis. In order to extract the structure functions from the 
equation for the amplitude, Eq.(^), we consider the components of the hadronic tensor. One 
can easily find from Eq.([|) that the structure function F 3 is determined by the ^-component 
of the hadronic tensor in the laboratory frame: 

7 F 3 = 2i W xy , (9) 

where by 7 we denote the ratio of the space to time componets of the momentum transfer, 

|q| / 4A/VV' 2 
7 * + — j ' <10) 

with x = Q 2 /2Mq . Using Eq.(§) we extract the nuclear F 3 from Eq.(^) and find, 

x A F 3 A (x A , Q 2 ) = -J- £ / [dp] x> Tr \A T {p) (^ {p, q) + £(p, g)/ T )l , (11) 

where xa = Q 2 /2Pa • q and x' = Q 2 /2p ■ q are the Bjorken variables for the nucleus and the 
bound (off-mass-shell) nucleon with four-momentum p, respectively. 

We stress here that Eq.(|TTD has been derived within a fully relativistic approach and 
no approximations have been made in terms of Q 2 . A similar approach (though focused at 
asymptotic Q 2 ) have been already applied to analyze the nuclear charged- lept on structure 
functions F lf2 |12[| as well as gi 2 |13[. One important observation which follows from Eq.flTTD 



is that, in general, traces do not factorize into completely separate nuclear and nucleon parts: 
the functions which correspond to different Dirac structures in Eq.([5|), enter F% with 
correspondingly different coefficients. We note also that these functions depend in general 
on three variables: x', Q 2 and p 2 . One arrives therefore to the conclusion that even in 
the impulse approximation nuclear structure functions are not determined by those of the 
on-mass-shell nucleon but are sensitive to their off-shell behavior. 

However in practice it may be quite sufficient to treat nuclei as nonrelativistic systems. 



The analysis of Refs. [T^JT^] shows that both for spin independent (Fi and F 2 ) as well as 



To be more precise A(p; Pa) is the imaginary part of the "hole" part of the full nucleon 
propagator. 
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for the spin (g\ and #2) nuclear structure functions one can in fact recover factorization in 
the nonrelativistic limit .0 In this limit nuclear structure functions factorize into convolution 
of nonrelativistic nuclear spectral function and "off-shell nucleon structure functions" . The 
latter are the definite combinations of off-shell inelastic formfactors (the functions in the 
discussed case) with a correct on-shell limit. 

In the next section we perform a similar nonrelativistic reduction of relativistic expres- 
sions in Eq. (O) to establish whether a similar factorization holds for the structure function 
F 3 . An important point of the present analysis is that it can be applied at arbitrary Q 2 . 



B. Nuclear Structure Functions in the Non-Relativistic Limit 

Our basic assumption in the remainder of the paper is that the nucleus is a non-relativistic 
system, made up of weakly bound nucleons. This necessarily involves neglecting antinucleon 
degrees of freedom, and corresponds to bound nucleons in the nucleus being slow, |p| <C 
M, \po — M\ <C M. In order to find nonrelativistic approximation to Eq.(]TT|) we have to 
perform a nonrelativistic reduction of all Dirac structures in the nucleon hadronic tensor. In 
doing so we follow the procedure which is described in detail in Ref. P"T|-|TB|]. Here we outline 
the basic steps of the method. The essential relation is the relation between the relativistic 
four- component nucleon field N and the nonrelativistic two-component operator ip (for the 
brevity of notations we suppress the isospin index r), 

where the nucleon operators are taken in a mixed (p, t) representation and the renormaliza- 
tion operator Z — 1 — p 2 /8M 2 guarantees baryon number conservation. We note here that 
Eq.(|12D is valid to order 1/M 2 in the series of 1/M-expansion and to this order does not 
depend explicitly on nuclear interactions (see discussion in Appendix of Ref. |I3]). We write 
the four-momentum of the bound nucleon asp = (M+e, p) and introduce the non-relativistic 
spectral function as follows: 

V(e,p) = J dte l£t (P A \^(v,0)i>{p,t)\P A ) I (P A \P A ) . (13) 

Notice that, as it follows from its definition, the spectral function is normalized to the 
number of protons (Z) and netrons (N) in the nucleus: 

J^V^(e, P ) = (Z,N). (14) 

By substituting Eq. flT2|) into Eq.(|^) we find the relation between relativistic and non- 
relativistic spectral functions and then reduce the four- dimensional Dirac basis to the two- 
dimensional spin matrices. We skip here technical details and summarize the results. The 



5 We note however the effect of mixing of the spin structure functions g± and 52 observed in Ref. 
& 
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non-relativistic limit of traces with the components of the operator Q^(p,q) can be written 
as follows: 



1 



2M A 
1 

2M A 



Tr 
Tr 



A(p)g (p,q)\ = V(s,p)F 3 (x',Q 2 ;p 2 ), 
A{p) g n (p,q)} = jjV{e,p) F 3 (x',Q 2 ;p 2 ) 



(15a) 
(15b) 



where n = x,y,z stands for the spatial components and F 3 (x', Q 2 ; p 2 ) is defined as follows: 

(16) 



F 3 ( X >, Q 2 ; p 2 ) = /W> ( 1 + 1 + ( 1 + ^=#^ + /« + /« 



2M 2 



M 2 



where for the brevety of notations we omit explicit dependence of functions /j on x', Q 2 and 
p 2 . Note also that F 3 in Eqs.flTo]) depends only on invariant variables. 

Few comments on Eqs. flTBD are in order. We emphasize here that non-relativistic limit 
is taken with respect to the nucleon momentum. In derivation of Eqs.(pT5|) we keep terms to 
order p 2 /M 2 and e/M and neglect terms of higher orders. At the same time Eqs.(15) are 
valid for arbitrary momentum transfer q. 

Eqs.(|15"D exhibit factorization of high-energy (structure function F 3 ) and low-energy 
(spectral function V) domains. We observe from ( |15|) that instead of four off-shell structure 
functions in non-relativistic domain we deal only with one specific combination of them, 
Eq.(16), which enters the nuclear structure function. We assume that Eq. flTB]) defines an 



analytical continuation of the nucleon structure function F 3 to the off-shell region in the 
vicinity of the mass-shell. It is easy to see that at p 2 = M 2 Eq . (|I6|) reduces to Eq.(||) that 
assures a correct on-shell limit of our analytic continuation of the structure function. 



C. Generalized Convolution At Asymptotic and Finite Q 2 



The factorization property of Eqs. ( |15"D allows us to write the nuclear structure function 
as a convolution of the off-shell nucleon structure function ([If]) and the nuclear spectral 
function (]T3l) . We substitute Eq.flTBD to Eq.flj]) and obtain the desired connection: 



x 



where we consider the nuclear xF 3 as a function of standard Bjorken variable x = Q 2 /2Mqo = 
xaMa/M and 7 is given by Eq . (flO|) . The variables which the off-shell nucleon structure 
function depends on read 

x 

X ' = l + (e + 7P,)/M' (18) 
p 2 = M 2 + 2Mb - p 2 . (19) 

We discuss now Eq. (|TT[) in more detail. We recall that no approximation is done with 
respect to Q 2 , so that Eq.(pT7|) is valid, in general, for any Q 2 . Let us consider first the limit 
of large Q 2 3> M 2 . In this case we put 7 = 1 in Eq . (|T7|) and find that xF 3 is smeared in 



S 



nuclei exactly like F2 fll2 |. However this similarity does not hold anymore when we study 
effects due to finite Q 2 



To explore Q 2 dependent nuclear effects in Eq.(jT7|) a useful observation is that all these 
effects come through the parameter 7. For the sake of convenience we consider the r.h.s. of 
Eq.flTTD as a function of 7 2 . Then expanding the under-integral expression in Eq.(|H]) about 
7 2 = 1 in Taylor series in r y 2 — 1 we obtain exactly the series in Q~ 2 . We skip here algebra 
details and show the resulting expansion up to Q~ 4 term, 



xF 3 A (x, Q 2 )/A «((1 + ^J x'F 3 (x', Q 2 ; p 

" UV %d* (x' 2 F 3 (x',Q 2 ;p 2 




M 

] X ' 1 G(x', Q 2 ;p 2 ) + %,x'd x ,G{x\ Q 2 ; p 2 ) ) . (20) 



Q2 I \ M V >r > M 2 



Here to simplify the notations the brackets (■ ■ •) denote the averaging over the nuclear 
spectral function, 



We also do not write explicitly the sum over the nucleon isospin states in Eqs.(20,2T[). Note 
that x' in Eq.(^) is given by Eq . (fl8|) with 7 = 1, and d x denotes the derivative with respect 
to x. The function G(x, Q 2 ;p 2 ) stays for the following combination of F3 and its derivative, 

G(x,Q 2 ;p 2 ) = 4xF 3 ( Xj Q 2 ;p 2 ) + x 2 d x F 3 (x,Q 2 ;p 2 ). (22) 



Few comments on Eq.(|20f) are in order. A natural question to ask is: What is an effective 
expansion parameter there? Formally this is M 2 x 2 /Q 2 . However this parameter is always 
multiplied by nuclear averages which involve powers of the nucleon momentum p z which is 
small in the nucleon mass scale. It may happen that this effect compensates a large factor 
M 2 . In fact it does. One can see this if one treats the Fermi motion and binding effects in 
Eq.(pOl) perturbatively, i.e. by expanding the structure function in (e +p z )/M about x' = x 
(detailed discussion of the expansion procedure see in Sec.[VI]). Then we find that the nuclear 
matrix element of Q~ 2 term is of order (p 2 )/M 2 . For the Q -4 term in the expansion Q2DD 
we find that two terms in the matrix element cancel each other to order p 2 /M 2 yielding the 
finite contribution only in order p 4 /M 4 . We conclude from this discussion that an effective 
expansion parameter in Eq.(pO|) is p 2 x 2 /Q 2 . 

We recall in this respect that Eq. (|l7D itself valid to order e/M ~ p 2 /M 2 . One must 



comment therefore that the nuclear Q 4 correction can not be calculated consistently without 
considering relativistic corrections to order p 4 /M 4 to Eq.([L7]). 



IV. OFF-SHELL NUCLEON STRUCTURE FUNCTIONS 

In order to obtain an estimate of off-shell effects in the nucleon structure functions we 
consider the latter in the leading order (impulse approximation for quarks, see Fig.^). In 
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analysing this diagram we consider the target momentum p to be generally off-mass-shell. 
As a convenient tool we use the Sudakov variables, i.e. parameterize the quark momentum 
in the loop, k, in terms of external momenta p and q, 

k = ap + f3q' + kj_, (23) 

where q 1 = q + £p and k± is two-dimensional vector (with k\ < 0) perpendicular to both p 
and q. We choose the quantity £ so that q' 2 = that gives 

^ = l + (l+4xV/Q 2 ) 1/2 ' ^ 

Note that £ is in fact the Nachtmann variable for the target with the mass squared p 2 . 

In doing the integration with respect to j3 we assume that the quark correlator has 
standard analytic properties in terms of the variables s = (k — p) 2 and k 2 , i.e. it is an 
analytic function in complex s- and fc 2 -planes with the cuts for s > and k 2 > 0. Using the 
dispersion relation for the quark correlator we can write the structure function as follows, 

F 3 (x, Q 2 - p 2 ) = J d ^ S f^ p(a, s, k 2 - p 2 ) 2p.qS ((k + qf) , (25) 

where < a < 1 and the integration over s runs over the spectrum of the "spectator" quark 
system and p describes the spectrum of these states.^ The squared quark momentum k 2 in 
Eq.(^5]) is 

k 2 = a(^-+p 2 ) + -^. (26) 
\a — 1 / 1 — a 

At asymptotic Q 2 ^> 1 GeV 2 one can neglect the quark virtuality k 2 (which is assumed 
to be finite) in the 5-function in Eq.(^) and the latter then goes to 5(a — x), so that we 
find in this limit: 

Fi s (x;p 2 ) = J^±p(x,s,k 2 ;p 2 ). (27) 

Now, we observe from Eq.(|27D that p 2 dependence of the structure functions has two pri- 
mary sources: explicit dependence of the quark spectral density on p 2 ("dynamical" p 2 
dependence) as well as p 2 dependence of the quark off-shellness k 2 ( "kinematical" p 2 depen- 
dence). Our general idea here is to fix the overall p 2 dependence of the structure functions 
by requiring the conservation of the normalization of F 3 off-shell. 

We recall in this respect that in the parton model F3 is interpreted as the valence quarks 
distribution in the target. The integral of F 3 , 

S G ls(Q 2 )= f'dxF^Q 2 ), (28) 
Jo 



6 To simplify notations we drop the explicit dependence of p on a the normalization point p 2 
which is fixed here to be Q 2 . 
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measures at asymptotic Q 2 the number of valence quarks (or the baryon charge) in the 
target, i.e. for the nucleon Sqls — 3 (the Gross-Llewellyn-Smith (GLS) sum rule f23f). The 
baryon current is conserved by strong interaction and, as a working hypothesis, we assume 
that S"gls is n °t renormalized off-shell (at least in the vicinity of p 2 = M 2 ). In the vicinity 
of p 2 = M 2 one can write this condition as follows: 

1 

^ dxd p 2F™{x;p 2 ) = 0, (29) 
where d p 2 denotes the derivative with respect to p 2 . In order to make use of this condition 



we will follow Ref. [12] and assume that the spectrum in s can be approximated by that 



calculated for a single effective mass s and write the spectral density as 

p = <5>(k 2 ,p 2 )5(s-s). (30) 

This model spectral density has been used in Ref. |12| to describe the valence quark dis- 
tribution in the proton. It was found that the k 2 dependence of (^0j) taken in the form 
(1 — k 2 /A 2 )~ n with the exponent n = 4 and the cutoff A 2 = 1.2 GeV 2 leads to a good 
description of the data at high Q 2 . To be specific we note here that these parameters have 
been fixed from the fit to the nucleon structure function data at average momentum transfer 
Q 2 « 5 — 10 GeV 2 , which fixes the normalization point for the spectral density. The leading 
Q 2 dependence of the structure functions can then be obtained by applying the evolution 
equation. 

Extrapolating p off-mass-shell we consider two possible models for p 2 -dependence in 
Eq. (|30"D . In the first model [13[] we assume that p 2 and k 2 dependencies factorize, $(/c 2 ,p 2 ) = 
ip(p 2 )§(k 2 ). Then the function (p(p 2 ) can be determined from Eq . (f29|) . In the second model 
|T2| we assume that the p 2 dependence comes from that of the cutoff, A = A(p 2 ). Its p 2 



dependence is determined again from Eq. (|29|) . 

Since we are interested in off-shell behavior of the structure functions in the vicinity of 
the mass-shell, the function of interest is d p 2F 3 taken at p 2 = M 2 . We found that both 
models give similar results for this function. Fig|| and Table | give an idea of the magnitude 
of the off-shell effect in the structure function F%. 

Going down to finite Q 2 we should comment that there is no direct reason to expect 
that Eq.(^) will remain valid, since S'gls does not correspond to the matrix element of a 
conserved current at finite Q 2 . Finite Q 2 corrections can be of two different kinds: dynamical 
ones, which are due to final state interaction of the struck quark with the spectator quark 
system; and kinematical corrections which are due to finite mass of the target. We focus here 
on the discussion of p 2 dependence of the structure functions. It is likely that the dynamical 
power corrections depend weakly on the target mass and we do not consider these effects 
here. In order to evaluate corrections due to finite mass of the target, let us consider the 
argument of the 5-function in Eq.(f25|) in more detail. Using Eq.(|23|) one can write it as 
follows, 

(k + q) 2 = 2p-q' ( a -0 + (k- &) 2 . (31) 

The last term in this equation measures the contribution from non-collinear components in 
Eq. (p3|) . Next we expand the 5-function in Eq.(|25|) in power series in (k — C,p) 2 - The zero 
order term in this expansion picks a = £ and we find the following contribution to F 3 , 
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(32) 



We found that the first order correction in (k — £p) 2 vanishes if the spectral density p does 
not depend on a explicitly (implicitly we have already used this assumption in Eq.fl3"0"|)). 
This can be seen after changing the integration d 2 k± — > dk 2 in Eq.(25[) and writing {k — £,p) 2 



in terms of the variables k 2 and s. We find therefore that the corrections due to non-collinear 
components in the 5-function in Eq.fl2~5|) appear starting from the order Q~ 4 and higher. It 
is then justified to use Eq.([3~2"l) to order Q~ 2 (though higher order terms could be important 
at x — > 1). To this end we expand Eq.(|32D in series in p 2 /Q 2 and keep only the first order 
correction. We then have, 

2 

12 2\ _ cias/„ . „2\ V „ Q / 2™/ 2' 



F 3 (x,Q 2 ;p 2 ) = F™(x;p 2 ) - ^xd x {x 2 F™{ X] p 2 )) (33) 

We find from the latter equation that to order Q~ 2 the right-hand side of Eq.(p3|) is 
expressed entirely in terms of the asymptotic part of F%. This observation allows us to 
calculate the derivative d p 2F 3 to order Q~ 2 using the asymptotic structure function given 
by EqsMM- 



V. NUCLEAR EFFECTS IN SUM RULES 



The moments of the structure functions are determined by the matrix elements of quark- 
gluon operators in the target. Those moments which correspond to the matrix elements of 
conserved currents are of special interest. We recall that F% obeys to the Gross-Llewellyn- 
Smith sum rule |23"|, i.e. the GLS integral Eq.(|28|) measures the baryon charge of the target. 



We stress however that this correspondence is approximate and exists only at asymptotic 
Q 2 . At finite Q 2 the GLS sum rule receives the QCD corrections. The following equation 
summarizes the theoretical status of perturbative as well as power twist-4 corrections to the 
GLS sum rule for the proton, 



3.25 



31.06 



a. 



27 Q 2 



M 2 f 1 
+ — - / dxx 2 F?(x,Q 
Q z Jo 



+ ■ 



(34) 



where a s = a s (Q 2 ) is the running QCD coupling constant. The coefficients of a 2 and ct^ 
have been calculated in |24| (the coefficients in Eq.(|4]) are given for Nf = 4). A power 
correction in Eq.(|3"4]) due to twist-4 operators has been evaluated in |2l)| , |2"6| by the QCD 
sum rules method. Ref. [25| gives ((O)) ~ 0.33 GeV 2 (with a typical QCD sum rules error 
about 30%). The estimate of Ref. |2§ is slightly larger, ((£>)) w 0.53 GeV 2 . As the last term 
in Eq.(|34D we write also explicitly the target mass correction. One observes from Eq.(p4|) 
a partial cancellation of the twist-4 power correction and the target mass term. We take 
((O)) of Ref. |26|] to obtain an upper estimate for the Q~ 2 term in Eq. (p4|) , then we have 
^GLs/3 = -0.11GeV 2 /g 2 . 
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Now we calculate nuclear correction to the GLS sum rule by integrating Eq.fl20D over x. 
In the right-hand side of Eq. (|20|) we change the integration variable x — > x' and after some 
algebra we have: 

S^) = (C*(l + ^ + i£M^±M>) ^ . (35 , 




We observe from Eq.(p5|) that nuclear corrections cancel out in the leading term in Eq.([20|) 
and we find this not surprising. Thus corrections appear starting from Q~ 2 order. We write 
explicitly also the Q~ 4 term just to illustrate our general conclusion about the expansion 
parameter in nuclear Q~ 2 series (see discussion at the end of Sec. [III C|) . This term introduces 



an extremely small correction and we do not consider it in what follows. 

Note that the overall nuclear Q~ 2 correction in Eq.(|35D comes from two sources: 1) the 
second term in the r.h.s. of Eq.(^) which is due to Fermi-motion of bound nucleons, and 
2) p 2 dependent Q~ 2 terms in the off-shell nucleon structure function, Eq . (|33"D . In order to 
evaluate the latter we observe that, since an average deviation from the mass shell in Eq.(RS]) 
is small, we can expand the off-shell nucleon structure function in Ap 2 = p 2 — M 2 about 
the mass shell point and keep only the first order correction. Using Eq.(|35|) and Eq.(p3|) we 
find, 

jS£ LS (Q 2 ) = S» LS (Q 2 ) + (Ap2 + |p2) f G dxx 2 F?{x,Q 2 ) 

+^ [ dxx 2 p 2 d p ,F 3 N (x, Q 2 ), (36) 

where we have assumed N = Z in the nucleus and F^ 1 = ^(F^+F 3 ) and Sq LS are respectively 
the structure function and the GLS integral for the on-mass-shell isoscalar nucleon. The 
derivative with respect to p 2 in the last term of Eq.(^) is taken at the mass shell point, 
p 2 = M 2 . The average off-mass-shell shift for the bound nucleon, (Ap 2 ) = (2Me — p 2 ), mean 
nucleon separation energy, (e), and mean nucleon momentum squared, (p 2 ), are determined 
by averaging with the nuclear spectral function by Eq . fl2~l|) . 

In order to evaluate the size of nuclear corrections in Eq. (|36"D , we refer here the results 
of computations of the spectral function for finite nuclei in the local density approximation 
f2"7| , |2"8"|j . In that approach the spectral function is decomposed in terms of a single particle 
contribution (mean field) and a correlation part. It is assumed that the latter is insensitive 
to the finite size of nuclei and is given by nuclear matter spectral function corrected with 
the help of local density approximation. As a result one finds that for a range of nuclei 
with 16 < A < 208 both the mean field and the correlation part of kinetic energy amount 



approximately 17MeV and slowly increase with A [28]. The mean field and the correlation 
parts of mean removal energies are about 30 MeV and 20 MeV respectively, and both increase 
with A. One of the traditional targets for neutrino DIS is iron 56 Fe. For this nucleus, basing 
on the discussion given in Ref. ||28|| , we take (p 2 )/2M = 35 MeV and (s) = —56 MeV. 
From these estimates we find (Ap 2 )p c = — 0.17GeV 2 . For comparision we give here also the 
corresponding values for the deuteron nucleus, (e)n = £d — (p 2 )d/2M, where en = —2.2 MeV 
is the deuteron binding energy. Taking the Bonn potential deuteron wave function we have 
(p 2 ) D /2M w 7 MeV and (Ap 2 ) D = —0.03 GeV 2 . Using these estimates as well as the results 
of Sec.|rV| we find from Eq.(RB) for the iron and deuteron nuclei, 
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5Sf,% 4.0- 10" 3 

<,LS (37a) 



3 Q 2 
5S5 TS 6.3 • 10~ 4 

-p = (37b) 

where the coefficients are taken in units of GeV 2 (5Sgls i n Eg. (|57|) denotes the nuclear 
correction to Sq LS in Eq.(|3(J), the factor 1/3 is introduced for the sake of convenience). 

It is useful to compare the sizes of different contributions to Eq.(|37]). The main con- 
tributions come from the second term in the r.h.s. of Eq.(|36|) which consists from the 
Fermi-motion term and the target mass correction term which are of opposite sign. Since 
the average off- mass-shell shift, (Ap 2 ), is negative and about as twice as large as the Fermi- 



motion term, (|p ), the overall contribution becomes negative. The last term in Eq . (pq) , 
which describes the deformation of the nucleon parton distributions due to off-mass-shell 
shift, brings about 20% of the total contribution to Eq.(^). Note that if we disregard the 
effect of p 2 -dependence of the bound nucleon structure function (which is a standard assump- 
tion in many convolution model calculations) we would leave only with the Fermi-motion 
correction which is about of the same size as qouted in Eq. (|37D but of opposite sign. 

We conclude from the present analysis that nuclear corrections to the GLS sum rule cancel 
out in the leading order in Q~ 2 series. A direct reason for this is that at asymptotic Q 2 the 
GLS sum rule corresponds to the baryon number of the target which is not renormalized by 
strong interaction. The GLS sum rule receives a tiny nuclear correction in Q~ 2 order which 
turns out to be more than an order of magnitude less than the corresponding QCD power 
correction. 

To conclude this section we give an expression for the arbitrary moment of nuclear 
structure function which follows from Eq. (|20|) . We integrate both sides of this equation with 
the weight x n_1 and keep terms to order p 2 /M 2 . We have, 

M n (Q 2 ) = f dxx n F 3 (x,Q 2 ), (38) 
jo 

+(A P 2 )d p iM^ (Q 2 ) + ^ (n + l)(n + 2)M» +2 (Q 2 ), (39) 

where d v iM n denotes the moments of d v iF 3 which can be calculated by Eq.(|33|) (see also 
Table H). We should comment however that application of Eq.(^) to large n should be 
taken with care. Indeed, Eq.(|3"9"D has been obtained by expanding the expressions like 
(1 + (e + p z )/M) n to order e/M ~ p\jM 2 . Obviously we are not allowed to do this when 
already the first order correction is of order of 1. Using the estimates for typical e and p 2 
from above we see that this happens for n ~ 10 for heavy nuclei. To calculate the moments 
with larger n one has to use the full nuclear spectral function. 



VI. DERIVATIVE EXPANSION FORMULA 



In this section we construct an approximation to Eq.(|i7|) taking further the advantage 
of a small parameter p/M. The smearing of the nucleon structure function in a nucleus 



14 



is governed by the nuclear spectral function V(e, p) which describes the distribution of 
bound nucleons over momenta and energies. Since characteristic momenta p and energies 
e are small as compared to the nucleon mass, an average deviation of the bound nucleon 
Bjorken variable, x', from the nuclear one, x, is also small. This simple observation allows 
us to obtain an approximate expression for the nuclear structure function in Eq.(|TT|) by 
expanding x'F^x 1 ', Q 2 ;p 2 ) in a Taylor series in p/M and e/M around p = (M, 0). In doing 
so we integrate the resulting expansion term by term and keep terms up to order e/M and 
p 2 /M 2 . As a result we obtain an approximation to Eq. (|T7|) in terms of a simple derivative 
expansion formula: 

xF 3 A (x, Q 2 )/A ~ xF 3 N (x, Q 2 ) - M xd x (xF 3 N (x, Q 2 )) + M x 2 d 2 x (xF 3 N (x, Q 2 )) 

+ <A P 2 >^ (xF 3 N ) (x,Q 2 ) + | M x 3 d?(x 2 F 3 N (x, Q 2 )). (40) 

Here we assume the isoscalar target with equal numbers of protons and neutrons and the 
coefficients are expressed in terms of averaged nucleon separation energy and momentum 
squared, Eq. ([21]) . 

We note that Eq. fljOf) being integrated over x reproduces Eq. (|36|) for the GLS sum rule 
as well as Eq.(|39|) for higher moments. In fact the inverse statement is also true, Eq.(f40|) 
can be derived from Eq . (|39"D by applying the Mellin transformation. 

The expansion formula (^) is similar to those which has been used in the discussion 



of the EMC ratio for the structure function F 2 (TT||T^j, as well as in the analysis of nuclear 
effects in the spin structure functions gi and g 2 [T3|] at asymptotic Q 2 . We note in this 
respect that ( ^0[ ) includes also Q~ 2 corrections due to nuclear effects. 

Figs.||| show the ratio R 3 = F^/AF^ calculated by Eq.(|40]) for the iron and deuterium 
nuclei for different Q 2 . Main features of the behavior of the ratio R3 calculated at large Q 2 
without the off-shell term are similar to those of the corresponding ratio for F 2 (EMC effect) 
and are determined by the competition of the separation energy and the kinetic energy 
terms in Eq.([E]). The position of the dip of R3 and its slope at large x are governed by the 
values of (e) and (p 2 ). The effect of the last term in Eq. fl4"0|) is similar to that of the kinetic 
energy term - it gives a positive correction. However this term is largely compensated by 
the off-shell term in Eq . (fiTf ) . We observe from Fig.f| a sizable correction due to the off-shell 



term. Its main effect is to increase the dip in Rs(x) and move its position to a bit larger 
x. Fig.[| shows the Q 2 dependence of the ratio -R3. Note that the main source of the Q 2 
dependence of the ratio is the off-shell term in Eq.(f40"D which incorporates also a target mass 
correction. We observe therefore that the size of the off-shell effects in the nuclear structure 
function even becomes larger the smaller Q 2 . 

We present also the estimates of nuclear corrections to the GLS integral truncated from 
below which is relevant in view of the extraction of the sum rule value from the data. 
Integrating Eq. (|40|) we have: 

Sgls(x,Q 2 )= C dyF 3 (y,Q 2 ), (41) 

J X 

S£ w (x, Q 2 )/A = S» LS (x, Q 2 ) + M x p 3 N ( x , Q 2 ) - M x %F 3 N (x, Q 2 ) 
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+ (Ap 2 )d p2 S% LS (x, Q 2 ) + l { -^r {-x 4 d x F 3 (x, Q 2 ) + 2 £ dyy 2 F 3 (y, Q 2 )) , (42) 

where c^Sqls denotes the GLS integral of d p 2F 3 . Fig|| shows the x dependence of nuclear 
corrections to the GLS integral. We observe that while the corrections practically cancel out 
at x = 0, according to the discussion of Sec.|V|, these are finite and negative in the region 
of practical interest and may reach up to ~ 10% for a heavy nucleus if one truncates the 
integral at x ~ 0.2. 

We conclude this section by discussing the region of applicability of the expansion formula 
(fl0|) which is both x and Q 2 dependent. The expansion procedure which leads to Eq.([40"D 
can be applied where the structure function is smooth. This is true when we are outside 
the elastic and resonance region, i.e. when (p + q) 2 > s with s ~ 2GeV 2 . The latter 
condition can be written in terms of the bound nucleon Bjorken variable as x' > x t h, where 
x t h = (l + (s — M 2 )/Q 2 )~ l describes the production threshold of the states with the invariant 
mass squared s. The bound nucleon Bjorken variable x' is smeared about the point x' = x 
by Eq.(|T%D. It is then clear that we do not get into a dangerous region if, 

1 - X/Xth ~ Pchar/M, (43) 

where p c har is a characteristic momentum of the bound nucleon. At asymptotic Q 2 , where 
x t h — ► 1, Eq.(|43"D tells us that the expansion formula Eq.(|40|) should be a good approximation 
for x 0.75 for heavy nuclei (for the deuteron the corresponding region is wider). We see 
also that for finite Q 2 the region of applicability of Eq.fl4"0|) becomes tighter. We note however 
that, as a matter of fact, even for large x outside the region (^) an approximate equation 
(|47]|) follows to the exact calculation by Eq.flTTI) though slightly overestimates the effect of 
Fermi-motion. This is not surprising because thresholds effects in the structure function 
restrict the integration region in Eq . (|i~7|) . 



VII. SUMMARY AND CONCLUSIONS 

In the present paper we have analyzed nuclear effects in the structure function F%. Our 
starting point was relativistic impulse approximation for the nuclear Compton amplitude. 

An important feature of the impulse approximation is the factorization property: nu- 
clear amplitude is factored into the nucleon Compton amplitude (generally off-shell) and 
a relativistic nucleon spectral function. However, because of off-shell effects, factorization 
between the amplitudes is not translated automatically into factorization between the struc- 
ture functions, that is assumed in usual convolution models. Indeed, as we discussed in 
the paper, the structure of the hadronic tensor of the off-shell nucleon is more complicated 
than that of the on-shell one. A direct reason for this is that one can not use the Dirac 
equation for the off-shell nucleon which greatly reduces the number of independent ampli- 
tudes in the hadronic tensor. In particular we found that the structure function F 3 splits 
into four independent inelastic form-factors for the nucleon off-shell. Similar observations 



have been already done for the structure functions F<i |12] as well as for g\ and #2 |[3 
Obviously these observations make more complicated the problem of extraction of nucleon 
structure functions from nuclear data and even could introduce a principal uncertainty into 
this problem. 
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We found that the factorization for the structure functions is remarkably recovered in the 
limit of weak nuclear binding when a nucleus is considered to be a nonrelativistic system. In 
order to observe this we have done a systematic 1/M-expansion of nuclear matrix elements 
which enter the nuclear Compton amplitude and found that to order 1/M 2 the nuclear 
structure function can be expressed in terms of the generalized convolution of the (non- 
relativistic) nuclear spectral function and a definite combination of the off-shell nucleon 
inelastic form factors which we take for the definition of the "off-shell nucleon structure 
function" . We note that this observation is generic for all structure functions p^2|JT3| . 

We stress here that our analysis is valid for arbitrary momentum transfer squared Q 2 . 
As a practical output of our approach we have developed a systematic expansion scheme of 
nuclear structure function in "nuclear twist" series, i.e. power Q~ 2 corrections caused by 
nuclear effects. We argued that an effective parameter in this expansion is p 2 har x 2 /Q 2 with 
Pchar being the characteristic momentum of the bound nucleon. This allows us to apply the 
Q~ 2 expansion down to small Q 2 ~ 0.1 GeV 2 and justifies to keep only Q~ 2 correction at 
Q 2 > 1 GeV 2 . 

One should comment however that these observations are done within the impulse ap- 
proximation and other sources of nuclear effects, such as the FSI, might be important in 
Q~ 2 order. It remains a challenging problem to quantify the FSI corrections in nuclear DIS. 

We applied our approach to calculate the moments of the nuclear structure function and 
in particular the Gross-Llewellyn-Smith sum rule. Nuclear corrections to the GLS sum rule 
cancel out in the leading order, which is due to the baryon charge conservation. We have 
calculated also nuclear Q~ 2 correction to the GLS sum rule which turned out to be negative 
and small (more than an order of magnitude less than the corresponding correction to the 
nucleon structure function itself as estimated by the QCD sum rules). Notice however that 
nuclear corrections are sizable for truncated GLS integrals that may have an impact on 
extraction of the value of the nucleon GLS sum rule from the data. 

We note here that in our analysis we can not avoid the problem of off-shell behavior 
of the amplitudes and have to deal with the problem explicitly. Unfortunately a rigorous 
approach, which would allow us to calculate the structure functions off-shell, is lacking and 
we have to relay more on model assumptions. In the present paper we evaluate off-shell 
effects basing on the non-perturbative parton model and assuming that the nucleon GLS 
integral is not renormalized off-shell, which is motivated by the baryon charge conservation. 
The latter condition allows us to evaluate the off-shell behavior of F3 at large Q 2 while the 
consideration of the target mass corrections leads to an estimate of this effect at finite Q 2 
as well. 

We found a sizable impact of off-shell effects on higher moments of the nuclear structure 
function and on its x dependence. Talking about the ratio of nuclear to nucleon structure 
functions R3 we observe that the off-shell effects tend to increase the dip at x ~ 0.6. This 
effect is even more prominent at finite Q 2 . In view of this observation we notice the impor- 
tance of the corresponding finite Q 2 analysis for the structure function F 2 as well as for the 
spin structure functions g\ and gi- 
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TABLES 



TABLE I. Moments of the off-shell term p 2 d p 2F^(x) calculated at p 2 = M 2 in the model 
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Sec. 
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P 2 d p2 M™ 











0.0394 


0.0165 
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p 2 d p2 M™ 




7.32 • 
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3.82 • 10" 4 
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1.11 • 10" 4 
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FIGURES 




FIG. 1. Dccp-inelastic scattering from a nucleus in the impulse approximation 




FIG. 2. Nucleon Compton amplitude in the leading order in Q 2 
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FIG. 3. The solid line shows the x dependence of the derivative p 2 d p 2(xF^)(x) calculated at 
the mass shell p 2 = M 2 in the model described in Sec. IV. Also shown is the nucleon structure 
function xF^(x) (dashed line). 
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FIG. 4. The ratio R 3 (x,Q 2 ) = \F^(x, Q 2 )/F^(x, Q 2 ) calculated for the iron (a) and the 
deuteron (b) nuclei by Eq.(|40|) at Q 2 = 5 GeV 2 . The dashed line shows the results without off-shell 
corrections to the structure function (i.e. without d p 2F 3 term), the solid line - full calculation 
including the off-shell term. 
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FIG. 5. Shown is the Q 2 dependence of the ratio R3 for the iron nucleus calculated by Eq. 
The dashed-dotted line corresponds to Q 2 = 15GeV 2 , the dashed line — Q 2 = 5GeV 2 , while the 
solid line — Q 2 = 3GeV 2 . 



1 . 1 



1 . 05 



. 95 



. 9 



. 85 



Sgls A (x,Q 2 )/A 
Sgls N (x, Q 2 ) 




FIG. 6. Shown is the x dependence of the nucleus/nucleon ratio of the GLS integrals, 
-^Sq LS (x,Q 2 )/ Sq LS (x,Q 2 ), calculated at Q 2 = 5GeV 2 . The dashed line corresponds to the 
deuteron, while the solid line to the iron. 
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